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ABSTRACT 



The X-ray source RXJ0648. 0-4418 is the only confirmed binary system in which a compact object, most likely a massive white 
dwarf, accretes from a hot subdwarf companion, the bright sdO star HD 49798. The X-ray emission from this system is characterized 
by two periodic modulations caused by an eclipse, at the orbital period of 1.55 d, and by the rotation of the compact object with a spin 
period of 13.2 s. In 201 1 we obtained six short XMM-Newton observations centered at orbital phase 0.75, in order to study the system 
during the eclipse, and spaced at increasingly long time intervals in order to obtain an accurate measure of the spin-period evolution 
through phase-connected timing. The duration of the eclipse ingress and egress, ~ 500 s, indicates the presence of an X-ray emitting 
region with dimensions of the order of a few 10* km, suiTounding the pulsar and probably due to scattering in the companion's wind. 
We derived an upper limit on the spin-period derivative |P| < 6 x 10^^^ s s^^, more than two orders of magnitude smaller than 
the previously available value. Significant X-ray emission is detected also during the 1.2 hours-long eclipse, with a luminosity of 
^3x10"^" erg s^^. The eclipse spectrum shows prominent emission lines of H- and He-like nitrogen, an overabundant element in 
HD 49798. These findings support the suggestion that the X-ray emission observed during the eclipse originates in HD 49798 and that 
the processes responsible for X-ray emission in the stellar winds of massive O stars are also at work in the much weaker winds of hot 
subdwarf s. 
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1. Introduction 

The X-ray source RX J0648.0-4418 is one of the very few 
known binaries composed of a hot subdwarf and a compact ob- 
ject. The optical component is the bright blue star HD 49798, 
classified as a subdwarf of spectral type 06, with effective tem- 

, per afiire Teg = 47, 500 K, a nd surface gravity logg = 4.25 ± 
0.2 ('Ku dritzki & Simonll978 ). Early spectroscopic observations 
showed a dominance of He and N lines and an underabundance 
of C and O, as well as radial velocity variations pointing to a 

• binary nature, that was later c onfirmed through the discovery of 

; an orbital period of 1.5477 d (lThackeravlll970b . 

The nature of the companion of HD 49798, undetectable in 
the optical band due to the presence of the much brighter sdO 
star, could be revealed only in 1996, thanks to the ROSAT dis- 
cover y of soft X-rays pulsed with a period of 13.2 s (llsrael et alJ 
Il997h . compatible only with either a neutron star or a white 
dwarf. While with the ROSAT data it was not possible to discrim- 
inate between these two possibilities, more recent observations 
carried out with XMM-Newton show ed that the compact ob ject 
is most likely a massive white dwarf ("Mereghetti et al."2009V 

For a distance of 650 pcQ ( Kudritzki & Simon 1978), the ob- 
served X-ray flux corresponds to a bolometiic X-ray luminosity 
of ~10"^^ erg s~^. Taking into account the rel atively well kno wn 
properties of the steUar wind of HD 49798 (lHamannll20Tol) . it 



can be shown that this luminosity is consistent with that expected 
for Bondi-Hoyle accretion onto a white dwarf, but too small for 
an accreting neutron star (iMereghetti et al.ll201 1) . The presence 
of an X-ray eclipse lasting '^1.2 h indicates that the system is 
seen at high inclination i ^ 82°. The masses of the two com- 
ponents have been dynamically measured through X-ray pulse 
timing and the spectroscopically determined optical mass func- 
tion: they are Mopt = 1-50 ± 0.05 Mq for t he subdwarf and 
Mx = 1.28 ± 0.05 Mq for its companion (IMereghetti et all 
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' The poorly constra ined Hipparcos parallax of 1.16±0.63 mas 
JPerrvman & ESAlll997h is consistent with this distance. 



120091) . Currently, HD 49798 is well within its Roche-lobe and 
mass transfer through wind accretion proceeds at a low rate of 
^8x10^^"^ Mq yr^ ^ . After the formation of a CO core, the sub- 
dwarf will expand and transfer He-rich material through Roche- 
lobe overflow at a higher rate. If the compact object is indeed a 

massive white dwarf, it may reach the Chandrasekhar limit and 

1 1 1 

either explode as a type la supernova (Wang & Han 2010) or 

form a fast-spinning neutron star through an accretion-induced 
collapse. This could be a way to create a millisecond pulsar di- 
rectly, i.e. without a spin-up phase in a low mass X-ray binary. 



During the only X-ray eclipse observed in May 2008 some 
faint X-ray emission was detected, but it could not be studied in 
detail with the short available exposure. We therefore caiTied out 
a series of time constrained XMM-Newton observations with the 
main objective to study the X-ray emission during the eclipse of 
this unique binary. 
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Table 1. Log of the 2011 XMM-Newton obsei-vations of RXJ0648.0^418. 



Observation 


Date 


MJD start 


MJD end 


Orbital Phase" 


Duration pn/MOS (ks) 


1 


2011 May 02 


55683.550 


55683.764 


0.685-0.824 


17.0/18.5 


2 


2011 Aug 18 


55791.876 


55792.069 


0.677-0.803 


15.0/16.6 


3 


2011 Aug 20 


55793.458 


55793.624 


0.699-0.808 


11.8/14.3 


4 


2011 Aug 25 


55798.040 


55798.267 


0.660-0.806 


18.0/19.6 


5 


2011 Sep 03 


55807.349 


55807.542 


0.675-0.801 


15.0/16.6 


6 


2011 Sep 08 


55811.997 


55812.189 


0.678-0.802 


15.0/16.6 



' Phase 0.75 corresponds to the middle of the eclipse. 



2. Observations and data reduction 

Six XMM-Newton observations of HD 49798/RX J0648 .0^4 1 8 
were performed between May and September 201 1 . Each obser- 
vation lasted about 15 ks and included orbital phase 0.75, cor- 
responding to the center of the X-ray eclipse (see Table [T] for 
details). Here we report on the results obtained with the EPIC 
(0.15-12 keV) and RGS (0.3-2.5 keV) instruments. The data 
were processed using version 1 1 of the XMM-Newton Standard 
Analysis Software. Some observations showed the presence of 
short time intervals affected by a high flux of soft protons, which 
caused a slight increase in the particle background at high ener- 
gies. We verified that our spectral results are insensitive to the 
inclusion or exclusion of such time intervals and therefore de- 
cided to use the whole data set. 

EPIC consists of two MPS and one pn CCD cameras 
dTurneretal .1120011: IStruderetalll200Ib . During all the obser- 
vations the three cameras were operated in Full Frame mode 
(time resolution of 73 ms and 2.6 s for pn and MOS, respec- 
tively) and with the medium optical blocking filter For the tim- 
ing analysis we used a circular extraction region with radius 
of 30". The photon arrival times were converted to the Solar 
System barycenter reference frame, by using the coordinates 
R.A. = 6'M8™04.7% Decl. = -44° 18' 58.4" (J2000), and 
corrected for the orbital motion of the source with the system pa- 
rameters given in Mereghetti et a l. (201 1 ). In order to avoid pos- 
sible contamination from a nearby source (located 70" away), 
we extracted the counts for the EPIC spectral analysis from a 
circle of 20" radius, which corresponds to an enclosed energy 
fraction of ^80% and ^75% for the pn and MOS, respectively. 
The spectra were rebinned so as to have at least 30 counts per 
energy channel. The background spectra were extracted from 
source-free regions on the same chip as the target. The spectral 
fits were done with version 12.7 of the XSPEC package. 



The RGS instrument (Iden Herder et al.lbOOlh has a signifi- 
cantly smaller effective area than EPIC, but its excellent spec- 
tral resolution makes it particularly sensitive to narrow emis- 
sion lines even in relatively faint X-ray sources like RX J0648.0- 
4418. We extracted the first order spectra and response matrices 
from both RGS units in each observation using a source extrac- 
tion region corresponding to 80% of the point spread function 
in the cross dispersion direction. This region is smaller than that 
typically adopted in RGS spectral analysis and was chosen in or- 
der to minimize the contamination from the nearby source men- 
tioned above (very likely an active star), which shows prominent 
X-ray emission lines. The source and background spectra and 
the response matrices of the seven observations were combined 
together to obtain a cumulative spectrum for each RGS camera. 
These spectra were rebinned in order to have at least 20 counts 
per energy channel. 




0.7 
Orbital Phase 

Fig.l. Light curve of RX J0648.0-4418 in the 0.15-10 keV energy 
range folded at the orbital period in bins of 535 s. The light curve has 
been obtained by summing the pn and MOS data of the six observations 
of 2011 and of the 2008 observation. The horizontal line indicates the 
background level. The vertical lines delimit the phase intervals used for 
the spectral analysis. 



3. Data analysis and results 

The 0.15-10 keV light curve obtained by folding the 2008 and 
2011 EPIC data (pn plus MOS) at the P = 1.547667 d orbital 
period is shown in FiglT] The horizontal line shows the level of 
the background and the vertical lines delimit the time intervals 
used for the spectral analysis described below. 



3. 1. Emission during tlie eclipse 

We start by examining the spectrum of the X-rays emitted dur- 
ing the eclipse, which can now be studied in detail thanks to 
the seven-fold increase in exposure time compared to the previ- 
ous observations. The spectrum of the single eclipse observed 
in May 2008 ('^4 ks exposure) was well fit by either a ther- 
mal bremsstrahlung (/cTbi— 0.55 keV) or a power law (pho- 
ton index = 2.8), but these simple models are rejected by the 
inclusion of the new data. In fact no simple single-component 
model can fit the total EPIC spectrum (2008 + 201 1, 0.2-lOkeV) 
extracted from the 0.735-0.765 phase interval (vertical dashed 
lines in FiglTli. In all cases the reduced x^ values are unaccept- 
able (2.8, 7.1 and 4.3 for absorbed power law, blackbody and 
bremsstrahlung, respectively) and the residuals from the best- 
fit model show the presence of a significant excess at 0.4-0.5 
keV, probably due to a blend of unresolved emission lines (lower 
panel of Fig EJ. A good fit is obtained with a power-law plus two 
naiTow line^ with energies fixed at Ei=0A3 keV and £^2=0.50 
keV, coiTesponding to emission from Nvi and Nvil ions. The 
best-fit spectrum of the pn data is shown in Fig|2] Fully con- 
sistent results were obtained with the MOS spectra. The best-fit 



" We fixed the width to cr = keV, which corresponds to the assump- 
tion of an intrinsic width smaller than the instrumental resolution. 
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Fig. 2. X-ray spectrum during the echpse (pn data). Top panel: data and 
best-fit model consisting of an absorbed power law and two Gaussian 
emission lines. Middle panel: residuals of the best-fit model. Bottom 
panel: residuals obtained using only an absorbed power law. 



parameters obtained by the joint pn and MOS analysis are given 
in Table |2] 

After background subtraction, the RGS spectra correspond- 
ing to the eclipse phases contain only about 100 counts in total, 
consistent with what expected from the best-fit EPIC model. An 
excess of counts is present at the energy of the most intense line. 
Fitting this excess with a naiTow line with energy fixed at 0.43 1 
keV, yields a line flux of (1.1 ± 0.7) x IQ-^ ph cm^^ g-i. No 
excess is seen at 0.5 keV, but this is not surprising considering 
the small photon counting statistics of the RGS spectra and the 
low flux estimated with EPIC for this line. 

To explore the possibility that the X-rays during the eclipse 
originate in HD 49798 we tried a fit with a combination of 
plasma emission models with different temperatures. The spec- 
tra of a large sample of O-type stars detected with XMM-Newton 
were well described by the sum of up to three thermal compo- 
nentfl with temperatures from ^0.1 to a few keV (Naze 2009). 
Using this model with plasma Solar abundances we could not 
obtain a good fit to the pn spectrum. This is not surprising 
considering the o verabundances of N and He in this system 
(iKudritzki & Simon 1978: Hamann 2010) . 

A good fit (reduced y^ = 0.86 for 23 d.o.f., Fig|3]l was in- 
stead obtained with the sum of three thermal plasma models with 
the He and N abundances set to the values observed in HD 49798 
(mass fractions of Xhc = 0.78 and X-^ = 0.025). The derived 
temperatures are kTi - 0.14, fcTa - 0.7, /cTg = 5 keXfl The 
0.2-10 keV flux is 7 x 10^^"' erg cm^^ s^^. This model is also 
consistent with the RGS data. 



3.2. Emission out of tlie eclipse 

We extracted from the six 201 1 observations the spectra exclud- 
ing the phase interval 0.725-0.775 (solid lines in FiglTJ. After 
checking that they did not show statistically significant differ- 
ences in flux or spectral parameters, we summed them to obtain 
a total spectrum for each camera. Assuming that the X-ray emis- 



Model MEKAL in XSPEC. 




The hottest component is required to account for the significant 
emission above E > 5 keV. Its temperature is very poorly constrained 
by the data, so we fixed it at the plausible value kTs — 5 keV. 



1 2 

Energy (keV) 

Fig. 3. Spectrum during the eclipse (pn data) fit with three thermal 
plasma emission models. Bottom panel: residuals from the best-fit 
model. 



sion observed during the eclipse is a stea dy component present at 
all orbital phases, as it was done also in lMereghetti et all (|20I Ih 
for the 2008 data, we include in the fit the eclipse model with 
parameters fixed at the values of Table |2] In addition, two other 
components are required to fit the data: a blackbody plus either 
a power law or a bremsstrahlung, yielding the results summa- 
rized in Table |3] Having included a more accurate description 
of the eclipse spectrum compared to our previous works, we ob- 
tain now slightly different best-fit values for the non-eclipsed 
emission. In particular , the blackbody component derived in 
iMereghetti et al.l (|20II|) had a higher temperature (^40 eV) and 
a smaller emitting radius (^18 km), compared to the updated 
values reported here. 

By summing the 2008 and 2011 RGS data corresponding 
to the uneclipsed emission we obtained a spectrum with total 
exposure time of 129 ks. This spectiTim suggests that the line at 
0.43 1 keV is present also out of the eclipse, since it shows an 
excess with flux of (0.65 ± 0.21) x 10^^ ph cm^-^ s^^ (i.e. a 
significance ~3(t). 

3.3. Eclipse duration and shape 

We fit the folded orbital light curve with a piecewise function 
consisting of two constants (the eclipse and out-of-eclipse val- 
ues) connected by linear functions to model the eclipse ingress 
(between ^\ and $2) and egress (between $5 and f&l). We as- 
sumed a symmetric profile, i.e. the same duration for the ingress 
and egress. With a fit to the pn light curve in the 0.15-0.5 keV 
energy range, where the signal-to-noise ratio is maximum, we 
obtained a total eclipse duration ($5 — $2) of 4311 ± 52 s and 
a duration of 525 ± 49 s for the eclipse ingress/egress. We re- 
peated the same analysis after dividing the 0.15-0.5 keV counts 
into two energy bands with roughly the same number of counts 
without finding statistically significant evidence for any energy 
dependence of the eclipse parameters. 

The gradual ingress and egress phases can be explained by 
either assuming that the X-rays originate from an extended re- 
gion around the compact object or that they are gradually ab- 
sorbed in the densest regions of the wind close to the sur- 
face of HD 49798. In the hypothesis of absorption in the atmo- 
sphere/wind of HD 49798, we would expect to observe some 
energy dependence in the orbital light curves near the eclipse. 
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As mentioned above, the comparison of the eclipse profiles in 
two energy bands was inconclusive, so we attempted to test this 
hypothesis by analyzing the spectra of the ingress and egress 
phases as follows. 

We summed all the EPIC/pn spectra extracted from the phase 
intervals 0.725-0.735 and 0.765-0.775 and tried to fit the result- 
ing spectrum by varying in two different ways the blackbody 
plus power-law model of the non-eclipsed emission. This was 
done a) by introducing a multiplicative factor (to simulate the 
gradual occultation of an extended source) and b) by keeping 
only the column density as free parameter (to check whether 
the flux decrease can be explained only by absorption). In both 
cases all the other parameters were kept fixed at their best-fit 
values of the non-eclipsed emission and the assumed constant 
eclipse spectrum (Table|2]i was included. The resulting values of 
reduced y^ = 0.99 and x^ = 1.42, respectively, favor the first 
interpretation. 

3.4. Spin period evolution 

The pulsations at the spin period of 13.2 s are clearly detected 
at consistent values in each of the six observations carried out 
in 2011. To improve the accuracy of the period determination, 
we performed a joint timing analysis of the 201 1 pn data in the 
soft X-ray band (<0.5 keV) as follows. We first phase-connected 
the two most closely spaced observations (n. 2 and 3) in or- 
der to obtain a preliminary solution. This was then iteratively 
improved by the successive inclusion of the other pointings ob- 
tained at increasingly longer time intervals (in the following or- 
der: 4, 5, 6, 1). The phase-connected ephemeris of the six ob- 
servations gave a period P2011 = 13.18424736(16) s, which 
is within the errors of the much less a ccurate value measured 
in May 2008 (P2008 = 13.18425(4) s. lMereghettiet al.ll2011l) . 
A linear fit to these two values plus all the other period mea- 
surements (2002 May and September with XMM-Newton and 
November 1992 with ROSAT) gives the following limits on 
the period derivative: -2.9 x 10"^^ < P < 3.6 x lO^^^ s s^^ 
(90% C.I.). We obtained a tighter constraint on P by extending 
the phase-connected timing analysis to the 2008 data and fit- 
ting the pulse phases with a quadratic function: the results are 
P = 13.1842474(2) s and |P| < 6 x lO^i^ s s^^. The 2011 
light curves folded at the spin period are shown in Fig. E]for 
the soft and hard energy ranges. As in the previous observations, 
the light curve at low energy shows a strong nearly sinusoidal 
modulation (pulsed fraction ^55%), while two peaks are visible 
above 0.5 keV. 



Table 2. EPIC spectral results for the eclipse emission. 



Table 3. EPIC spectral results for the uneclipsed emission. 



Parameter 


Value 


ATh (lO"" cm- 2) 


<2.9 


Photon index 


1 07+0.22 


P^L(10""ergcm-2s-i) 


7.3 ±0.7 


£1 (IceV) 


0.43 (fixed) 


h (IQ-^piicm-^ s"^) 


1 1+0.7 


E2 (IceV) 


0.5 (fixed) 


/2(10-^plicm-2s-i ) 


0.68«;^^ 


d.o.f. 


35 


xl 


1.16 







BB+PL 


BB+Brem 


iVn (10^" cm-2) 




1.6 ±0.7 


0.5 ±0.5 


fcTBB (eV) 




30.6 ±1.5 


33.5 ±1.3 


Rbb (km) 




6et% 


35tf2 


Photon index 




l-96triS 


- 


fpL(10-^-'*ergcm-2s- 


-') 


i.oit^:;]^3 


- 


fcTsr (keV) 




- 


3.7tg:? 


FBr(10-"ergcm-2s- 


') 


- 


O.lStVot 


d.o.f. 




262 


262 


2 
Xv 




1.27 


1.42 



" Blackbody emission radius for d = 650 pc. 
'' Flux in the range 0.2-10 keV. 




Phase 

Fig. 4. EPIC pn light curve of RX J0648.0-4418 folded at the spin pe- 
riod in the soft (0.15-0.5 keV, top panel) and hard (0.5-10 keV, bottom 
panel) energy ranges. The dashed lines indicate the background level. 



4. Discussion and conciusions 

The presence of X-ray emission also during the eclipse, already 
noticed in May 2008, is confirmed by the new observations of 
HD 49798/RX J0648.0-4418, which show a spectrum character- 
ized by strong N emission lines. A possible explanation is that 
these X-rays originate in HD 49798 itself, i.e. they are unrelated 
to the presence of a compact object in this system. This is par- 
ticularly interesting since no 0-type subdwarf has been detected 
in X-rays up to now0. 

Our current understanding of the X-ray emission from early 
type stars is based on OB main sequence, giant and supergiant 
stars, which, compared to HD 49798, have very different lumi- 
nosity, mass and composition. The X-ray emission of these stars 
originates in their strong, radiation-dri ven winds, where plasma 
is heated by shocks and instabilities dPallavicinil Il989l) . High 
resolution spectra of bright O type stars are characterized by 



" Flux of the power-law component in the 0.2-10 keV range. 



^ with the exception of BD+37°442 JLa Palombara et all l20T2h . 
However, the possible presence of pulsations at 19 s suggests that also 
BD +37° 442 has a compact companion. 
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the p resence of numerous emission lines (see, e.g. iKahn et aP 
I2OOI ). but data with lower resolution and/or statistical quality 
are usually well fit by a combination of plasma emission mod- 
els with different temperatures. The temperatures we derived 
are fully consistent with those fo und in a large sample of O 
stars studied with XMM-Newton (lNazel l2009). The luminosity 
of '--^3 X lO''" erg s^^ that we measured during the eclipse corre- 
sponds to an X-ra y-to-bolome tric flux ratio of ^10^^, as typical 
O-type stars ( Pall avicinil 1 989h . Thus we conclude that the X-ray 
luminosity and spectrum observed during the eclipse are consis- 
tent with those of early type stars and support the interpretation 
of HD 49798 as the first hot subdwarf detected at X-ray energies. 

In this respect it is interesting to note that HD 49798 is one of 
the few s dO stars for which e vidence of wind mass loss has been 
reported (iHamann et al.lfl981l) . This suggests that the processes 
responsible for X-ray emission in the stellar winds of massive 
O stars are also at work in the much weaker winds of hot sub- 
dwarfs and makes HD 49798 a particularly interesting target to 
explore the properties of radiatively driven winds through X-ray 
observations over a wide range of stellar parameters. 

Thanks to the better coverage of the eclipse orbital phases, 
the new observations allowed us to study the eclipse shape. We 
derived a duration of 43 11 ±52 s for the total eclipse, consis- 
tent with the rough value of 4300 s that was estimated from 
the 20 08 observation and used to constrain the system inclina- 
tion in' Mereghettietal.1 (l2009h . We could also estabHsh that the 
eclipse ingress and egress occur gradually. This can be explained 
if part of the X-ray emission originates in an extended region 
surrounding the compact object. It is easy to show with a sim- 
ple geometrical model, based on th e known orbital parameter s 
and radius of HD 49798 (1.45 Rq, iKudritzki & Simonlll978h . 
that the observed duration of ~500 s of the echpse ingress and 
egress implies dimensions of the order of a few lO'' km, much 
larger than the radius of a massive white dwarf (^3000 km). 
The presence of strong pulsations, which must necessarily orig- 
inate close to the compact object, sets a limit on the flux coming 
from the extended region. Considering that ~55% of the counts 
below 0.5 keV are pulsed, this limit is of about 2 x 10~^^ erg 
cm~^ s~^ (0.2-0.5 keV). A plausible scenario is that part of the 
high-energy radiation emitted by the compact object is scattered 
in the wind of HD 49798, as observed in more luminous X-ray 
binaries p owered by wind accretion, such as, e. g.. Vela X-1 or 
Cen X-3 dSako et al.lll999t lEbisawa et al.lll996l) . This interest- 
ing possibility deserves to be explored with more sensitive X- 
ray observations with good spectral resolution (Sartore et al. in 
preparation). 

The new observations confirm the stability in the long term 
properties of the X-ray emission from this binary system: the 
luminosity and pulse profile do not show significant variations 
compared to previous observations, but thanks to a better model- 
ing of the eclipse spectrum we could obtain more accurate values 
for the spectral parameters of the uneclipsed emission. The lower 
temperature of the blackbody-like component, which dominates 
the low-energy pulsed emission, implies an emitting area diffi- 
cult to reconcile with origin on the surface of a neutron star. 

Thanks to a properly planned spacing of the observations, 
we could obtain for the first time a phase-connected timing so- 
lution for this system. This provides stronger constraints on the 
secular evolution of the 13 s spin period: the new limit on the 
period derivative, |P| < 6 x 10^^^ s s^^ (90% c.l.), is more 
than two orders of magnitude smaller than the previous value 
(-5x1 0"" s s^i < P < 9x10"" s s"\ iMereghetti et all 
(1201 ih ). Considering the small torque expected for the low ac- 



cretion rate onto this pulsar, this limit is still insufficient to defi- 
nitely rule out a neutron star, a possibility not dismissed yet, al- 
though we believe that the observed X-ray luminosity and spec- 
trum are best explained by a white dwarf companion. In this re- 
spect it is essential to extend the phase-connected solution over 
a longer time span by means of new X-ray measurements of the 
spin period. 
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